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Abstract. Motivated by the interest in models of the early universe where statistical 
isotropy is broken and can be revealed in cosmological observations, we consider an 
SU{2) theory of gauge interactions in a single scalar field inflationary scenario. We 
calculate the trispectrum of curvature perturbations, as a natural follow up to a previ- 
ous paper of ours, where we studied the bispectrum in the same kind of models. The 
choice of a non-Abelian set-up turns out to be very convenient: on one hand, gauge 
boson self-interactions can be very interesting being responsible for extra non-trivial 
terms (naturally absent in the Abelian case) appearing in the cosmological correlation 
functions; on the other hand, its results can be easily reduced to the U{1) case. As 
expected from the presence of the vector bosons, preferred spatial directions arise and 
the trispectrum reveals anisotropic signatures. We evaluate its amplitude t^vl, which 
receives contributions both from scalar and vector fields, and verify that, in a large 
subset of its parameter space, the latter contributions can be larger than the former. 
We carry out a shape analysis of the trispectrum; in particular we discuss, with some 
examples, how the anisotropy parameters appearing in the analytic expression of the 
trispectrum can modulate its profile and we show that the amplitude of the anisotropic 
part of the trispectrum can be of the same order of magnitude as the isotropic part. 
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1. Introduction 

The Cosmic Microwave Background (CMB) radiation is a window on the early universe 
physics [1] and represents a major object of the WMAP [2, 3, 4] and Planck mission [5] 
searches and of other future cosmological investigations. 

A powerful feature of the CMB is its degree of non-Gaussianity [6, 7]. A signal is con- 
sidered Gaussian if the information it carries is completely encoded in the two-point 
correlation function, all higher order connected correlators being zero. Present mea- 
surements have already allowed to put interesting limits on the level of non-Gaussianity 
of the CMB, and the recently launched Planck satellite and future experiments are 
expected to provide either tighter constraints or possibly a detection of a primordial 
non-Gaussian signal. Many of the existing inflationary models are in excellent agree- 
ment with the CMB observations but differ from each other in terms of predicted amount 
of non-Gaussianity [8, 9, 10, 11, 12, 16, 17, 19, 20, 21, 22]. If we want to be prepared for 
an accurate comparison with the data in the near future, it is therefore crucial to under- 
stand what the different models forecast in terms of non-Gaussianity. This translates 
into calculating three, four and higher-order correlation functions for the cosmological 
fluctuations. 

In this paper, we focus on the trispectrum of curvature perturbations [23, 24, 25]. 
This has been calculated in several models, such as multifield slow-roll inflation 
[26, 27, 28, 29], in the curvaton mechanism [30], in theories with non-canonical ki- 
netic terms both in single field [31, 32, 33, 34] and in the multifield case [36, 37]. % 
We consider here another class of models, specifically those including primordial vector 
fields [13, 41, 42, 43, 44, 45, 46, 47, 48, 49]. The interest in these models finds its mo- 
tivations, besides in the search for a strong non-Gaussian signature, also in the study 
of statistical anisotropy. In fact, several anomalies in the CMB data suggest signatures 
of statistical anisotropy breaking [50, 51, 52, 53, 54]. This unexpected feature could, 
among other ways, be explained with the existence of preferred spatial directions de- 
fined by some primordial vector fields, in models where they can somehow leave their 
imprints on the cosmological fluctuations. 

Recent works where two and three point functions of the curvature perturbation are 
computed in the presence of primordial vector fields include Refs. [46, 47, 48], which 
consider a U{\) gauge theory and different specific models (vector inflation, vector cur- 
vaton, time-dependent gauge coupling models); Ref. [45], which works in the Abelian 
case and for a model of hybrid inflation; a previous paper of us, [55] , where an SU (2) 
gauge theory is explored. The interest in extending the existing work to a non- Abelian 
case, was for us mostly driven by the possibility of new contributions to the correlation 

X At present there arc no constraints on the non-linearity parameters characterizing the trispectrum 
from CMB analyses [38, 39]. The only constraint available on the parameter (7nl of "local" cubic non- 
linearities has been obtained from Large-Scale-Strucuture (LSS) data yielding I^nlI < 10^, and future 
CMB and LSS observations could probe I^nlI > 10"* [40]. 
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functions arising from self-interactions of the gauge fields in the Lagrangian. These new 
contributions are naturally absent in the Abelian case. For the bispectrum, we proved 
that they are potentially comparable and, in some case, even larger than the Abelian 
contributions. The same motivation leads us to study, in this paper, the trispectrum for 
the same kind of theories. In fact one of the reasons why the study of the trispectrum 
is important is that in most cases the three and four-point correlation functions are 
correlated to each other, being underlined by the same coupling constant. 

Throughout the paper, we will employ the SN formula [56, 57, 58, 59], in order 
to express the fluctuations in the e-folding number in terms of the fluctuations of all 
of the quantum fields on the initial temporal slice. In addition to that, the Schwinger- 
Keldysh formalism [60, 61, 62] will be used for calculating the non-Abelian contributions. 

The organization of the paper goes as follows: in Section 2 we introduce and briefly 
discuss the Lagrangian of the model; in Section 3 we list all the terms contributing 
to the trispectrum of the curvature perturbation, classifying them in three separate 
categories, scalar, vector and mixed (depending whether they originate from scalar, 
vector or both fields); in Section 4 we calculate the mixed terms; in Section 5 we are 
faced with the vector terms, particularly with the non-Abelian contributions from the 
gauge field trispectrum, represented by vector-exchange and point-interaction diagrams; 
in Section 6 we provide an estimate of the t^l parameter; in Section 7 we study the 
shape of the trispectrum; in Section 8 we discuss how the anisotropic features modulate 
the trispectrum and show that the amplitude of the anisotropic part can be of the 
same order of magnitude as the isotropic one; in Section 9 we extend our calculations 
to models of gauge interactions where the kinetic term is multiplied by a function of 
time; in Section 10 we present our conclusions. The Appendices provide more details 
about some lengthy calculations and expressions shortened in the rest of the paper for 
the mixed contributions (Appendix A), for vector-exchange (Appendix B) and point- 
interation diagrams (Appendix C). 



2. The model 



Several models have been proposed which include primordial vector fields. We want 
to consider a scenario of a scalar-field driven inflation, where the inflaton field coexists 
with an SU{2) gauge multiplet. The latter can play a non-negligible role in determining 
the amount of curvature perturbations through mechanisms that will be later discussed 
in the paper. A possible choice for the Lagrangian is 



S 
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a=l,2,3 



^ 1 a=l,2,3 



(1) 
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being the scalar field Lagrangian and F^^, = d^B^ — dyB'^ + gcS^'^'^B'^^B^. The physical 
fields are defined as = (^-Bq, B"'/a{t)^ and ^ is a numerical factor. The kinetic term 
for the gauge fields is multiplied by a generic function / that can be viewed as a function 
of time. 

The gauge field quantum fiuctuations are expanded as follows 



SAnS^v)= fi%e^'-' E K(g)<^''M^(g,r/) + ef(-g)(a-J)^Mr(g,r/)], (2) 

l^^J X=L,R,l 



where a and are creation and annihilation operators and the index A runs over left, 
right and longitudinal polarizations. 

Depending on the functional form of /, on the value of the bare mass mo and of C,, 
different models arise. In the special case / = 1, it was proven in [13] that, for a small 
gauge field mass mo ^ H and choosing ^ = 1/6, the transverse components of the 
gauge field perturbations satisfy equations of motion similar to the one of a minimally 
coupled light scalar field which, therefore, has a an almost scale-invariant spectrum on 
large scales. In this case, the gauge fields acquire an effective mass = mg — 2if^ (it is 
correct to use the approximation R ~ —12H^ during inflation). The reader can refer to 
Appendix A of [55] for all the details of the equations of motion for the background and 
for the gauge field perturbations and their solutions. On the other hand, the longitudinal 
wave function in this model has been shown to behave like a ghost, its kinetic energy 
being negative, and it is still debated whether or not it has a well-defined quantum field 
theory, also because it is governed by an equation of motion suffering from singularities. 
In spite of this, a regular solution was proposed in [46]. The physical validity of this 
model was questioned in [65, 66, 67, 68]. However these instabilities do not represent 
an issue in some models with a varying kinetic function where / (and eventually the 
effective mass) is instead not a constant [45, 63, 64]. 

In view of these different possibilities, we will employ for the longitudinal mode the 
parametrization in terms of the transverse mode 55" = n{x)5B'^ we introduced in [55], 
motivated by the intent of keeping our calculations as general as possible, in spite of 
having to make a specific choice for the Lagrangian. Our calculations can indeed be 
easily specified for different models. We are going to first consider the Lagrangian with 
mo ^ H and ^ = 1/6, keeping / = 1. In the last section of the paper, we show how to 
extend our calculations to /(t) models. 

3. Trispectrum from 5N formula 

We want to calculate the trispectrum of the curvature perturbation 

{CuSkSkSk,) = i^^fsHh + h + h + hmh, h, h) (3) 

using the 6N expansion. In the presence of primordial vector fields [46, 55] 



C(x , t) = N^6<P + N^6A^^ + ^-N^^ {6<pf + ^A^JM^M^ + N^M^A^ 



(4) 



where is the number of e- foldings between an initial time t*, which is often 
conveniently fixed at the horizon crossing era of the given wave number, and the final 
time t, labeling the slice of observation for (. The partial derivatives of have been 
thus defined 



'dN\ 



(5) 



\d<pdAl)^, 
and so on for higher order derivatives. 

Once Eq.(4) is plugged in Eq. (3), several terms will arise. It is convenient to collect 
them in three categories: (purely) scalar (S), mixed (M) and vector (V) 



(6) 

We anticipate here that T^'^ and TY can be written as the sum of an Abelian, T^'X, 



and a non-Abelian, T, 



■c 

M,V 
((NA) 



contributions. This will help stress what the origin of the 



different terms is and help derive the non-Abelian limit of our results. 

Calculations of the purely scalar contributions already exist in the literature, 
whereas the mixed and the vector terms have not been computed so far. We report 
here the tree-level result [27, 28] 



NfT^{ki,k2,h,k4) 



+ N^N^^ [P^{ki)B^{\ki + /csl, /cs, k^)+ perms. 
+ NlNl^ [P^{k,)P^{k2)P^{\k, + perms. 
+ iV|iV000 [P4>{ki)P^{k2)P4,{k^) + perms] . 



We defined 



(7) 



= (27r)3<5(3)(fci + h)Pm^ (8) 

(^hMkMk) = (27r)'5(=^)(^i + k2 + h)B^{ki, k2, A^s), (9) 

(%,%/0a;3%4) = (27r)35(=^)(^i + k2 + h + h)T^{ku 4, h), (10) 

where, in single-field slow-roll inflation, to leading order we have: P^ = H'^/{2k^) {H^ 
being the Hubble parameter evalutated at horizon exit k = aH); from [12, 9] we learn 
the result 

VeH^M{k,) 



ir 



mp 



M being a function of the momenta moduli ki of dimension (mass) ^; finally from 
[26, 29] we have 

T^c^HlM{h), (12) 
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where M is a function of dimension (mass) 

The new contributions from vectors will be evaluated in the next sections. 



4. Calculation of mixed contributions 

We list the total (connected) contribution from the mixed (scalar-vector) terms 



C ~ 4)^^ (l>a^^ <i>b 



P^tih + h)PJh)PJk2)+ perms. 



+ N^,NI:NINI, PZ{h)P::t{k2)PA\h + perms 



P^{h)P^{h)P^;,{h)+ perms. 



lip 



perms. 



+ N^^N^N^N^ P^{k2)P^{\h + k2\)P;,:{k^)+ perms. 

Pa'c'ih)Pb7iki + k2)P4h) + perms. 

+ NINI^N;, [p^l{h)B^{k,, k2, \h + h\)+ perms.' 

+ N^^Nl^N^N^^ [P^iks)Bf:^ih, h, h + h) + perms. 
+ [higher order terms], 

where we defined as usual 



(13) 



(14) 

{SAI,:^6AI,:^6A;,:J = {27rf6^^\h + h + h)B;%ih, h, h)- (15) 

In Eq. (13), we avoided including terms that involve mixed scalar- vector bispectra, the 
reason being that in the chosen Lagrangian no direct coupling exists between scalar 
(infiaton) and vector degrees of freedom. Before proceeding with the calculations, one 
more remark should be made which will hold for the rest of the paper: the partial 
derivatives of the e-folding number w.r.t. the temporal components of the gauge fields 
will be set to zero, the reason being that the are equal to zero (see Appendix A of 
[55]). 

The power spectrum of the gauge fields was previously calculated in [55] and it reads 

(16) 



ptm 



ab 
long 



where 



(k) 



rpeven 
ij 

Tlp%k) 



^ijkkk 

kj^kj . 



ki kj , 



and 



P^" =5ab6A''^ik,t*)6A''*{k,t*), 
Pf =6abSAUk,t*)6Ai*{k,t*), 

Plong = ^abSA1^j^g{k,t )6 A''i^j^g{k, t 
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(17) 
(18) 
(19) 

(20) 
(21) 
(22) 



= {1/2){P^ ± Pl^), SA'^ j^i^^g indicating the eigenfunctions for the gauge fields. 
The total power spectrum of ( then becomes [55] 

P^{k) = P''\k) [l + g'^'^ (k ■ Na) {k ■ N,) + is'^'k ■ (Na X iV,)] , (23) 

where g'^' = {P^tg - Pf)l[NlP^ + (iV, ■ N^) Pf] and s^' = P^l'/lN^P^ + (iV, ■ iV,) P^^] 
and the isotropic part of the power spectrum is 



P'-^'^ik) = NlP^{k) + (n, ■ N,) Pl^ (24) 

As to the vector modes bispectrum (in the last line of Eq. (13)), this was also introduced 
and computed in [55] using the Schwinger-Keldysh formalism (see Eq. (50) therein) and 
it turns out to be 



h) = glHl e'^ac^abe^e f f dx) T^.T^.T^ + pevms. + c.c. (25) 

\ K J ^^^^^ \J / a/37 

We define x = —kr], where k = ki + k2 + k^ + k^. As a reminder, (/ dx)^i^^ indicates the 
time integral over the wavefunctions of the internal legs multiplied by the wavefunctions 
of the external legs of the diagrams, the indices a, /3 and 7 runnig over longitudinal 
and transverse polarization states only. In the model we adopted, no violation of parity 
occurs, therefore P'^ = and P^* = P^. 

The next step is to calculate the time integrals in Eq (25). It is well-known that, for 
the Lagrangian (1) with / = 1 and an effective mass = ttlq — 2H^ ^ —2H^ the 
transverse mode behaves as a light scalar field during infiation [13] 



SB' 



(26) 



2v^ 

Throughout the paper, we find it convenient to adopt the parametrization introduced 
in [55] for the wavefunctions of the longitudinal mode in terms of the transverse mode, 
i.e. 

6B\\ =n{x)6B^ , (27) 

where n is un unknown function of a; = —kr]. Like in [55], the choice of this 
parametrization is motivated by the convenience of a more general approach when 
handling longitudinal modes, in view of their still debated nature and of the possibility 
of making our analysis adaptable to more than one physical model. 
When performing the time integrals, we should remind ourselves that, for the kind 
of wavefunctions we are dealing with, the major contribution comes from the region 
of integration around the horizon H ^ k/a [9]. Also, let us assume that n{x) is a 
smooth function: this is reasonable, based on requirement that the longitudinal mode 
asymptotically scales as the transverse mode at late times so as to produce a scale 
invariant spectrum. For example this happens in the models with the transverse mode 
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given by (26) with mo — [13] , or in models with varying kinetic function and mass 
[63, 64]. § It then appears to be a good approximation evaluating n{x) at the horizon 
crossing time x* and taking it out of the integrals in Eq. (25), which leads to 

f dx) = - ^,-,^-121212 *^ l^^EEE + {Beee COS X* + Ceee sin X*) EiX*] , 
J J EEE 24:k-^k(k2k^x*^ 

dx) = [ dx) = [ dx] = n^(x*) i dx) , 

/EEl \J / ElE \J / lEE \J J EEE 

dx] = i dx) = i dx) = n'^{x*) i dx) , 

/HE \J / lEl \J J Ell \J J EEE 



dx) = n^{x*) { dx) . (28) 

/ III \J J EEE 

The letter 'E' is an abbreviation for 'even', whereas '1' stands for 'longitudinal' and the 
definition of the functions Aeee, Beee and Ceee can be found in Appendix C of [55]. 
The final expression for line 8 of Eq. (13) becomes 

T^'lunes = N,P^{k,)glHl (^)' (/ dx)^^^ [iPee + n\x*) [i^ei + 41 + I^'i^e) 

+ n\x*) (/fi + Jg! + Jgi) + nV)/S^] (29) 

where the expressions of the anisotropy coefficients r^f^^ can be found in Appendix A 
(Eqs. (A.l) through (A.4)). This is the only non-Abelian contributions to T^^, the re- 
maining terms are also present in the Abelian limit Qc ^ ^■ 

In some vector field scenarios, the mixed scalar-vector derivatives A^^^' vanish, so 
Eq. (13) does not contribute to the trispectrum. This can certainly happen in those 
models where a direct coupling between scalar and vector fields is missing, but the 
latter condition is not sufficient for concluding that the mixed derivatives are null. As 
an example, it is useful to refer to [69] , which, among other things, includes an analytic 
study for the case of a set of slowly rolling fields with a separable quadratic potential. 
The number of e-folding is written as a sum of integrals over the different fields, to be 
evaluated between their values at an initial (generally set at around horizon crossing) 
and a final times. For each field, the value at the final time depends on the total 
field configuration at the initial time, so the mixed derivative of A^ can in principle be 
non-zero. Anyway, if the final time approaches the end of infiation, it is reasonable 
to assume that, by then, the fields have stabilized to their equilibrium value and no 
longer carry the memory of their evolution. If this happens, the sum of integrals which 
defines A^ becomes independent of the final field configuration and its mixed derivatives 
can therefore be shown to be zero. It turns out that we are allowed the same kind of 
analytic study, if we work with the Lagrangian in Eq. (1) and if we introduce some 
slow-roll assumptions for the vector fields (see [55] for details, in particular Appendix 

§ The assumption that ~ is crucial if by "late time" we mean a epoch within the end of inflation 
and well-after horizon crossing, since in the opposite case, i.e. for non-negligible mg, the solution (26) 
for the transverse mode would no longer apply and it would become harder to produce a scale-invariant 
power spectrum. 
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A for a discussion about what these assumptions imply and Appendix B for the actual 
calculation of and its derivatives). 



5. Calculation of purely vector contributions 

Let us now turn to the vector trispectrum contributions: they are quite similar to the 
scalar ones (Eq. (7)) 



= KN;:NpNST^t;iih, h, h, h) 

+ Ni^Nl^N^N^^ [pf^{h)B':^{k + 4, ^3, k,) + perms'. 
+ N^N^,N^^^,Nl] [P^^;{h)Ptt{k,)P^f^{k^ + h)+ perms 



perms. 



(30) 



where 



(31) 



The power spectra and the bispectrum appearing in Eq. (30) were provided in the 
previous section (Eqs.(16) and (25)). All we are left with is then calculating the 
trispectrum of the gauge fields, i.e. T^^^^ in the first line of Eq. (30). However, in 
the next pages, we are also going to rewrite the other lines of the expression above for 
completeness and for the computation of the non-Gaussianity parameter tnl in Sec. 6. 
Like in [55], the terms that are proportional to Nq will not be taken into account since 
A'^ = for all a = 1, 2, 3. 

5.1. Ahelian terms in 

They include the third and fourth lines of Eq.(30). Let us rewrite them in such a way 
that their anisotropy coefficients are easily readable. 

Let us first introduce the following quantity (we use the same notation as in [55]) 



M',{k) = KPtk'ik) = Pl\k) [5,uK+p''\k)h (k ■ Na) + iq'^^k) {k X Na 



(32) 



where p°-'^{k) 
Also, let us define 



JDac 
long 



pac\ /(^pac^^ qac(j^^^ = pac^pac g^^^ jy^ ^ {N^,N^,N^). 



Li{k) ^ N^:,Ptl{k)N'^ 



= PfimL ■ Nif+p'f{k) [k . {k . Ni^) + tq'^{k)k . Ni, X Nif] 
where iVi,^(A^i], AT;?], AT;?^^). 

Using the equations above, the third and the fourth lines of (30) become 



rpV 



M^L^MJ + MiMJM^Nf, 



(33) 



(34) 



where the subscript A always stands for "Abelian" and A^^gd third derivative of 

A^ w.r.t. the vector fields. Notice that the anisotropy coefficients p"^ and g"'^ become 
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zero respectively for theories in which transverse and longitudinal modes evolve exactly 
in the same way and P"* = 0. 

In the situation where the vectors Na are all aligned along a unique spatial direction, 
some more considerations can be made about the level of statistical anisotropy in 
Eq. (34). In particular, for theories where q"''' = applies, 7"^^^ isotropizes if the 
four wave vectors kt lie in a plane perpendicular to the direction of the gauge vectors; 
otherwise, for theories where p"''' = 0, ^^^a) ^i^^ appear isotropic when considering wave 
vectors parallel to the gauge vectors. 

5.2. Non-Ahelian terms in 

The non-Abelian contributions to T'^ are given by the first and second line of Eq.(30). 
The bispectrum -B^^^ was reported in Eq.(25). If we define A'^^ = MfN^^ (where Mf 
was introduced in Eq. (32)), then the second line of (30) can be rewritten as 

Triune 2 = KNlN^B^^^ih + h) (35) 

which is formally similar to Eq.(50) of [55], with one of the replaced by N^. The 
result is therefore given by 

T^liine 2 = glH^ (^^) dx^ [Ieee + n^(x*) (^Ieei + hiE + Iiee) 

+ n\x*) (luE + iiEi + hi) + n\x*)im] + perms. (36) 

where the lap-y are provided in Appendix A after Eqs. (A.1)-(A.4). It is interesting 
to notice that, in those models where the longitudinal and the transverse modes have 
the same time-evolution (i.e. n{x) = 1) and parity is preserved (i.e. P"'' = 0), the 
contribution from Eq. (36) becomes isotropic 



Triune 2 = V'^V^'^AT, . N^N, ■ A' (37) 



k 



X 



24k^klkik'lx*^ 



{Aeee + {Beee cosx* + Ceee sinx*) EiX*) + perms 



Another observation worth to be made is that the /q,/3^ become all equal to zero in the 
case where the three gauge vectors point in the same spatial direction. In order to 
prove this, one should employ, besides their definitions in (A.1)-(A.4), also the explicit 
expressions for the derivatives of A^ (which can be found in [46] or also in [55] and that 
will be used in Section 6 in the evaluation of tnl)- 



We are now left with one last term to be analysed, i.e. the first line of Eq.(30). 
Like for the other non-Abelian contributions, this will be calculated using the Schwinger- 
Keldysh formula 



(fi|e(t)|fi) 







T 



jJlHi{t')dt' 



e,(t) 



T e 



-if^Hi{t')dt' 







(38) 



10 



Figure 1. Diagrammatic representations of vcctor-enchange (on tlie left) 
and contact-interaction (on the right) contributions to the trispectrum. 



where \Q) and |0) are the vacua respectively of the interaction and of the free theories. 
Both time-ordering, T, and anti-time-ordering operators, T, are needed to account for 
the fact that the quantity on the left-hand side is not a scattering amplitude between 
an initial and a final state, it is instead an expectation value of a given cosmological 
observable at a given time. The subscript / reminds the reader that we are working in 
interaction picture, so all the fields can be treated as free and thus expanded in terms 
of creation and annihilation operators in the usual way. 

To leading order in the SU (2) coupling and in the perturbative expansion of the gauge 
fields, there are two types of diagrams, the vector-exchange and the point interaction 
diagrams (see Fig. 1), which arise respectively from third and fourth-order interactions 
in the Hamiltonian 



Hlnl = 9ce''"'9''9''d,SB]6Bl6B^ (39) 
H^^l = gy^'e^^''g^^g''5B^5Bl5B'^.5Bf. (40) 

Both of these contributions to TY are expected to be of order ~ qIH^. 



5.2.1. Vector- exchange diagrams 

Let us begin with the two vertex diagram. Using the language of Eq. (38), it can 
be put in the form 



dv iH^v)-H^iv ) ) (41) 



A '- J—oo ^ 

where now H = H^^t^ © = 6Alj^6Al6Ap6A'^ and the inclusion symbol as usual points 
out that what stands on the right-hand side is only one of the contributions to (0(?7*)). 
The superscripts -|- and — refer to different rules of field-operator contractions, i.e. 



6Bt'^ivyB'^'^iv") = t^tHv'^vMv' -v") + f^tHv\v"Mv" -v' 

6Bt''-ivlSB';'-iv") = Tli;iv\v"). 
6BrivlSB'^'\v") = t^tHv^v). 
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In Fourier space we have 



ng'(^) = T^^'''{k)P+ + iT°f\k)P^^ + Tl;;^%k)Pf^ (42) 



Iif^{k) = T^p''{k)Pl' + tT^f''{k)Pt^ + Tlp'{k)Ptl' (43) 

where = {1 /2){P^ ± P^^) , P'^ being equal to the product of the two eigenfunctions 
6abSBf{k,r])6B*j^\k,r]') and Pf = (similar definitions hold for Pf and Pif^g). 

We remind the reader that in our model Tff'^ = 0. 

Eq. (41) can be rewritten as follows 

{Qiv*))^^^{T[Q{A + B+C+D)]) (44) 



where 



A = r dr]'H+{r]') P rfr/"i7+(r/") 

J—oo J—oo 

B= r d7]'H-{r]') r dr]"H-{7]") 

J—oo J—oo 

C=- r dr]'H+{r]') r dr]"H-{7]" 

J—oo J—oo 

D=- r dr]'H~{r]') T dr]"H+{7]" 

J—oo J—oo 



For each one of the integrals listed above, due to the presence of both the fields and their 
spatial derivatives in H^^^, there are three different sets of contractions of the external 
with the vertex field-operators: for the first set, the field-operators with derivatives in 
the vertices are both contracted with external fields; for the second one, only one of the 
two field-operators with derivatives contracts with an external field (the other contracts 
with another internal field); for the third set, the field-operators with derivatives contract 
with each other. 

A sample set of contractions of the first type is provided in the following equation 

u J 'J ' "J' " 

rpabcd —^ Jc b c a b c i i 



X 



I dr^'a\ri) j d^\\ri')g^^g^ ^' g^'^g^ (45) 

where the first four lis correspond to contractions between external and internal fields 
whereas the last one indicates the contraction between the two remaining internal 
fields. The factor comes from expressing the external (physical) fields in terms 
of the comoving ones. As a reminder, we define {2Ti:)^5^^\ki -|- ^2 + ^3 + ^'dTij^f = 
{5A15A^^5A1.5A^). The expression in (45) can be rewritten as follows 



rpabcd 9c 



I I I II n II I I II II f II 

_^abc^a b c k^k^,S'"' 5^^ S"" 5'^'' " 



X 



12 



in jm im nn 



(46) 



where, again, the greek indices of the sum indicate either the transverse (E) or 
longitudinal (1) modes, the J drj' J drj"^ stand for the integrals over the wave functions, 

the chosen time variable being x = —kr] {k = Z]j=i,...,4 h = \ki\). 
Let us define the coefficients T°^t7^" = kmk'T"Tl^^T?',Tf,T'',. They should be 
calculated for each one of the three different sets of contractions and for each permutation 
within the specific set. This is a straghtforward but rather lengthy and not particularly 
interesting calculation. A convenient way to proceed could be the following: we first 
compute the time integrals in order to find out which one among the combinations 
of longitudinal and transverse mode functions in the string [a, /3, 7, 5, a] provides the 
highest amplitude for the trispectrum (in order to be able to perform this comparison 
we work, as it is usually done when trying to quantify the amplitude of a three or 
of a four-point function, in the so called "equilateral configuration", which for the 
trispectrum means taking k\ = k2 = kj, = kii = fc^g = ^14); for the combination with the 
highest amplitude, we then calculate the coefficients T°!j^'i^" for all the different sets of 
contractions and sum over all the permutations. || 

Let us now perform our calculations. The wavefunctions we are going to adopt were 
introduced in Sec. 4 (see Eqs. (26) and (27)). It is possible to verify that B = A* and 
D = (T and that integrals of type A are consistently smaller in amplitude than integrals 
of type C. We therefore report the combined contribution C + D = 2Re[C\ for one of 
the permutations 



dt] / d7] 



EEEEE 8klk^kiklk^-^{kf2 + h + ^2) (^12 + k3 + k^jx 



*8 



X [(M - 2E)[{N - 2F){AB + CD) + {2H + L){CB - AD)] 

+ {2G + P)[{N - 2F){AD - CB) + {2H + L){AB + CD)]] (47) 



dr] 
drj 
drj 
drj 
drj 



drj 
drj 
drj 
drj 
drj 



EEEEl 



EEEll 



EElll 



EUU 



UUl 



2i * 

n [x 



4 / * 

n [x 



n^ix* 



n^(x*) 



drj 
drj 
drj 
drj 
drj 



drj 
drj 
drj 
drj 
drj 



EEEEE 



EEEEE 



EEEEE 



EEEEE 



(48) 
(49) 
(50) 
(51) 
(52) 



EEEEE 

where A, B, C, D, E, F, G, H, L, M, N and P are functions of x* and of the momenta 
moduli to be provided in Appendix B (Eqs. (B.l) through (B.9)). Obviously the value of 
the integrals does not change when permuting its labels a(5'j6a, besides having a different 

II Notice that this procedure provides an approximate result, but it is certainly allowed since, when 
working in the equilateral configuration (as we do when we evaluate the amplitude of the trispectrum) , 
it is easy to realize that the integrals (^J drj J drj ^ are independent of the specific set of operator 
contractions, which only differ by the T^^^"^^"^ and by the Levi-Civita coefficients. 
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power of the coefficient n{x*) multiplying them. We need now to find out if there is 
one, among the integrals in Eqs. (47) through (52), that has the largest amplitude, 
i.e. understand if something can be said about the order of magnitude of n{x*). We 
could try to extrapolate some information about n{x*) from what happens at very late 
times. In the models discussed in Ref. [13, 46] it turns out that the longitudinal mode is 
= \/25B^.^lf this is the correct asymptotic behaviour and we find it reasonable to 
extrapolate back until the horizon crossing epoch, it is then correct to conclude that in 
this case the amplitude is the largest for the integral among the ones listed in Eqs. (47) 
through (52) containing the highest powers of n, i.e. for J drj' J ^f]"^ The coefficients 
we intend to calculate are then of the kind T/jj^" only. We list them below for the three 
different sets of contractions we introduced above and for one particular permutation 



^ kihiki2M {ki ■ kf^^ (ki ■ (ks ■ h) ■ kf2) > (53) 
rpim^) ^ kskf^kusA {ki ■ (^2 ■ ^2) (^3 ■ ^4) (^3 ■ ^12) , (54) 
rj^iiim ^ k^^k^^ki2u {ki ■ ^2) (^2 ■ kf^ {h ■ ki2) (^4 ■ k^^^ ■ (55) 

We adopted the following notation: ks = \ks\^ kg = ks/ks, the index s running over 
the four external momenta; kss's"s"' = ksiks'jks"kks"'i, with s,s',s",s"' = 1,2,3,4 and 
with the indices k, I indicating the spatial components of the vectors; fc^-, = ks + kg', 
k^y = \ks + ks'\ and so k^^, = fc/^,//c/^,. 

It is possible to prove that, once the Levi-Civita coefficients and the sum over the 
permutations are taken into account, only the first set of contractions is left (see again 
Appendix A for detailed calculations). The final result after these cancellations can be 
written in the following form 

{SA^SA'jSAlSAf), D {2nf6^^\h + ^2 + 4 + h)gl (^)\^'''"e'^''" \l x ^1234 x (j^U 



u=l 



12 



k 

+ 11 X ki324 X + /// X kii32 X U j J . (56) 

All the possible permutations have been included in the previous equation and, as a 
reminder, the indices k, I are hidden in kss's"s"' on the right-hand side. We define 



*8 



8kfkiklklkl^^{kf2 + ki + k2){k^2 + h + h)x 
X [(M - 2E) [{N - 2F) {AB + CD) + {2H + L) {CB - AD)] 

+ (2G + P) [{N - 2F) {AD - CB) + {2H + L) {AB + CD)]] (57) 

(from Eqs. (47)- (52)). The function // is defined from / by exchanging /c2 with k^ and 
ki2 with k^2,'i similarly, /// is defined from I by exchanging k2 with k/^ and kf2 with k^^, 

^ As another example, in models with varying kinetic function and mass [63, 64], we have verified that 
n{x) ^ 1 at late times and for a vector field that is light until the end of infiation. 
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so they are all functions of the horizon crossing time x* = —krj* and of the moduli of 
the external momenta and of their sums. This amounts to seven independent variables, 
X*, fci, ^2, ^3, ^4, and /cf4. The coefficients U {i = 1, 12) come from T^^?^'' and 
so they are also functions of the momenta moduli (see Eqs. (B.IO) through (B.21) for 
their expressions). Finally, the anisotropic part of Eq. (56) is represented by the kss's"s"' 
terms, which, in the final expression for the curvature perturbation trispectrum, have 
their spatial indices contracted with the derivatives A^^" of the number of e-foldings w.r.t. 
the vector fields as follows 



(C^kSk^Ck^CkJ ^ N^N^mNf{6At6A';6Al6Af). D i27rf6^'\h + h + h + h)9. 

X I X ti^ X A/ + perms. 
where the anisotropic term in the first permutation is 



* s 4 



It can be interesting to rewrite Aj in terms of all its variables 

, 2 



(5^ 



(59) 



3 

{a<b)a,b=l 



(60) 



:i,j] = [l,2],[3,4] 

The M/'s are 2x2 matrices whose entries are represented by the cosines of the angles 
between the wavevectors and the N°- 



cos 6ia COS 6ja 



COS 9ib COS 6 



'jb 



i.e. cos^ia = ki ■ N"- and so on. 

The two permutations in Eq. (58) can be written in a similar fashion with anisotropic 
coefficients 



[i,j] = [l,3],[2,4] 

{N^f[N'yx n det[Miif'') . 

[ij] = [l,4],[3,2] 

The number of angular variables is equal to 12. These are to be added to the six 
scalar variables from the isotropic part of (58) {ki, /c2, k^, k^^, k^2 ^^^d k^^ and to three 
parameters represented by the lengths of the vectors N"- in A. 

Like in (36), the anisotropy coefficients here become equal to zero in the event of an 
alignment of the gauge vectors along a unique direction. 

Notice that the isotropization we found out to occur in the n{x) = 1 case for the 
trispectrum contribution in Eq. (36) (and that will be as well proven for the contact 



A,, = e"''' e""' N^N^jNlNtkn2i = J2 

(a<b)a,b=l 



Am = e''''' N^N!;N'^Nfku32 = E 

(a<b)a,b=l 
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interaction in the next subsection), does not hold for the vector-exchange diagram. 
In fact, even without having to calculate the T^^l^" coefficients explicitly, it is easy 
to verify that, summing all of them up, we are left with a term proportional to 
^ahc^ade (^g^ . . jy^^ fj^^ . j\f^^ ^ comlug from T^Jf^^^, plus its permutations. The 

anisotropy therefore survives because of the presence of spatial derivatives in the 
interaction Hamiltonian that are contracted with the external field operators. 

5.2.2. Point-interaction diagrams 

Let us now move to the one-vertex diagrams 

(e(r/*))Dz(T[e r (61) 

where now H = h\^^ and again = 5A'^^5A\,5A'p5A'^. After working out the Wick 
contractions, this becomes 

xY.([dx) T^.T^^T^J'^.+permutatzons. (62) 



" mnop mi nj oi pj 



Let US list the coefficients = '^mi^nj'^oi^rri ^'^^ permutations 



+ A;3434 + k3232 — ^^2 ' ^4 (^3234 + ^1214) ~ ^np^ml^ol + ^1414 + ^1212 

+ Snpkmikosh -ks-ki-ks (A;i232 + A;i434) + h ■ ■ hi, (63) 



'^mnt^^ — ^mokpAkni — ^mo^n2^p4^2 ' ^4 ~ ^3434 + ■ (/i;3234 + ^1214) ~ ^1414 

+ A;i434A;i ■ h - ki234ki ■ hh ■ h, (64) 

'^mnop — ^npkoskmS — ^upkmlkoskl ■ k'i — A;3434 + /Ci ■ /C3 (A;i434 + /C1232) — ^3232 

+ ^3234^:2 ■ ki - ki234h " hk2 ■ h, (65) 

"^mnop — ^mokn2kp2 — 5mokn2kpik2 " ^4 — ^3232 + ^^2 " ^4 (^3234 + ^1214) — ^1212 

+ ki232ki ■ h - A;i234A;i ■ A;3A;2 ■ k^, (66) 

"^mnop — kmlkolSnp — kmlkoaklsSnp — /Cl414 + ^1 " ^3 (^1434 + ^1232) " ^1212 

+ kuuh ■ ki - ki234ki ■ hh ■ h, (67) 

'^mnop — ^3434 ~ ^3234^2 ' ^4 ~ ^1434^1 ' ^3 + ^1234^1 ' ^3^2 ' ^4; (68) 

^mnop 5mokn2kp4:k2 ' ^4 k2 ■ ki (^3234 + A;i2i4) + A;i234A;i ■ A;3^2 ■ ki, (69) 

'^rnnop — ^3232 — ^3234^2 " ^4 " ^1232^1 " ^3 + ^1234^1 " ^3^2 " ^4, (70) 

'^mnop — ^1212 — ^1214^2 " ^4 — ^1232^1 " ^3 + ^1234^1 " ^3^2 " ^4, (71) 

'^mnop ~ ^1414 ~ ^1214^2 ' ^4 ~ A;i434fci ■ /C3 + /i;i234^1 ' ^3^2 ' ^4; (72) 

^ranop ^npkmlkoskl ks- ki- k^ (A;i434 A;i232) + /i;i234/i;i ■ ^3^2 ■ ^4, (73) 

^™ = ^1232^1 ■ A;3 - ^1234^)1 " ^3^:2 " L (74) 



'^mnop ~ ^1214^2 ' ^4 ~ ^1234^1 ' ^3^2 ' ^4; (75) 



" mnop 

nlEU 



"^mnop — ^1434^1 ' ^3 ~ ^1234^1 ' ^3^2 ' ^4; (76) 



- mnop 
-lElU 



^mnap = ^Z2Z^h2 ■ - k^2?,^k^ ■ k^k^ ■ k^, (77) 

T^mnop = kusiki ■ kh ■ h. (78) 

When evaluating the integrals (/ dx)^i^^g, we use again the wavefunctions previously 
introduced in Eqs. (27) and (26). The final result is 

/ dx] = ^,,.,2,\,2,2 *7 [Qeeee + Aeeee cix* {Beeee cos X* + Ceeee sin x*) 

J J eeee 2Ak^kik2kikix*' 

+ Deeee six* {Eeeee cosx* + Feeee sinx*) ], (79) 
dx] =n^{x*)( f dx) , (80) 

/ EEEI \J / EEEE 

dx) = n\x*) (fdx) , (81) 

/ EEll \J ) EEEE 



dx = n^(x*) { dx] , (82) 

/ lllE \J / EEEE 

dx) = n\x*) (fdx) , (83) 

/ V"' / EEEE 

where Qeeee, Aeeee, Beeee, Ceeee, Deeee, Eeeee and Feeee are functions of x* 
and of the momenta ki = \ki\, ci and si stand respectively for the Coslntegral and the 
Sinlntegral functions. The expressions of these functions can be found in Appendix C. 
It is again important noticing that the anisotropy coefficients become zero if the gauge 
fields are all aligned. 

Finally, summing up the coefficients in Eqs. (63) through (63), one realizes that if the 
longitudinal and the transverse mode evolve in the same way, the total contribution 
from the point-interaction diagram is isotropic 



^ o^^5^2^2^2^2 ^j [Qeeee + AEEEEcix* {Beeee cosx* + Ceeee sinx*) 
+ Deeees^x* {Eeeee cos x* + Feeee sin x*) ] + permutations. (84) 



6. Calculation of t^l 

The non-Gaussianity of a given theory of infiation and cosmological perturbations can 
be studied by looking at the expression of the two well-known parameters fjqi and Tjqi. 
The parameter rjvi, is defined from the trispectrum normalized with the isotropic part 
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of the ( power spectrum + 

^ 2T(^{ki,k2,h,h) 

pi.o(^^)piso(^2)P^^°(A;fJ + 23 perms.' ^ ^ 

where P**° was provided in Eq. (24), and we recall that k^^^ = ki + k^^. 
The total tnl will then in general be a function the vectors ki, i = 1, ...,4 and can be 
written as the sum of scalar and vector contributions (we remind the reader that we are 
working under the conditions that the mixed contributions are equal to zero) 

r^L = rS + rffi. (86) 

More explicitly, we have 

(5) _ 5,1 , 5,2 , 5,3 , 5,4 /g^s 
'^NL — '^NL + '^NL + '^NL + '^NL^ \° ' ) 

(V) _ V,l , V,2 , V,3 , VA /ggN, 

where we counted, both for the scalar and the vector parts, four different contributions 
which can be read from lines 1 through 4 of the right-hand sides of both Eqs. (7) 
and (30). Notice that and t^I are the non-Abelian contributions, t^'^ and r^'f^ 
the Abelian ones. More precisely, r^'^^ refers to the contribution from the vector field 
bispectrum, Eq. (36), refers to the vector- field trispectrum and therefore includes 
both vector-exchange and contact-interaction terms. Let us now carry out a quick 
estimate of the different contributions to tnl, ignoring vector and gauge indices. The 
scalar contributions are 

rfjl - (89) 

5,2 _ 5,3 _ 5,4 _ € /Qp,\ 
"^NL — '^NL — '^NL ~ p^3- ) 

where we define f3 = {Na/Nip)'^. The vector contributions are included in the Table 1 
for the general case (without specifying the explicit expressions for the derivatives A^^ 
and Naa), in vector inflation and in the vector curvaton model. 



Table 1. Order of magnitude of the vector contributions to tnl 
in different scenarios. 



V,2 



V,3 

'''nl 



V,4 

'''nl 



general case 



in3 f 'gc 



10 



9c 



o3/2 3/2 



V. inflation 



same as above 



10- 







V. curvaton 



same as above 



10 (1+^)3 [h, 



2ILE.V 

A ) 



[LE.V 
A ) 



The previous expressions were derived using ~ {mp^y 



m 



p ' 



N, 



p, with e 



^)/ (2mpif^). Numerical coefficients of order one in the amplitudes 



+ This normalization choice, which simplifies the expressions of tjvl compared to normalizing with the 
complete from Eq. (23), is motivated by the fact that the anisotropic part of the power spectrum 
cannot noticeably affect our results, being weighted by parameters such as g"^^, that have to be kept 
much smaller than one in order to avoid an anisotropic expansion. 
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were not reported. The 10^ factor in r^'^ is entirely due to the vector-exchange 
contribution, which turns out to have an amphtude that is several orders of magnitude 
larger than the one of the point-interaction term. For the derivatives of w.r.t. the 
vector fields, the results of [46, 55] were employed, which we schematically (neglecting 
all indices) report below 



— , Naa ^ — , 
nip rrip 



Naaa = (91) 



in vector inflation, and 



Na^j, Naa^^, Naaa^^ (92) 

for the vector curvaton model, where r = (Sp^) (Sp^ + 4p0) (p^ and indicating re- 
spectively the energy densities of vector fields and inflaton at the epoch of the curvaton 
decay). Here what we call "vector inflation" is a regular scalar-field- driven inflation oc- 
curring in the presence of vector fields whose contribution to the total universe energy 
density is assumed to be subdominant w.r.t. the inflaton energy density in such a way 
that the expansion remains approximately (and within experimental bounds) isotropic. 
The vector curvaton model, proposed in [13] for the first time, is similar to the classic 
curvaton mechanism [14, 15, 16, 17, 18, 19], except for vector fields, rather than scalars, 
playing the role of the curvaton field. 

The parameters (3, e, Qc and r appearing in the previous expressions are all as- 
sumed to be smaller than unity. Concerning /3, it is well-known in the U{1) case that it 
has to be kept small in order to prevent too much anisotropy from affecting the power 
spectrum (see [54] where data analysis is provided in this sense; this is also discussed 
in [46]); no data analysis has up to now been performed considering more than one 
preferred spatial direction, so, unless all the gauge vectors are about aligned with one 
another, the upper bound on /3 is not known. It sounds nevertheless safe and reasonable 
to assume that this quantity is quite small also for a generic orientation of the vector 
fields in our SU (2) model. The parameters e and gc are obviously supposed to be small 
respectively for ensuring slow-roll of the inflaton field and allowing perturbation theory 
to hold. Finally, in the curvaton model, r has to remain small at least until the end of 
inflation for an almost isotropic exponential expansion to occur. 

The ratio mp/H is of order 10^ (assuming e ~ 10~^ and from the observed power spec- 
trum of fluctuations); as to the other ratios A/H and rrip/A appearing in Table 1, their 
value is not strictly constrained. In fact, if on one hand it is reasonable to assume that 
the expectation value of the background gauge fields does not exceed the Planck mass, 
in principle no constraint can be put on A/H, which can generically be a very large 
number, even though specific gauge field configurations exist where it is of order one. 
An example of such configurations is the one where all the gauge fields have the same 
magnitude and are not approximately aligned (see Section 6 and Appendix A of [55] for 
a complete discussion). To conclude our observations on Table 1, we can say that r^'l 
and Tjy'^ can be either smaller or much larger than one both in vector inflation and in 
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the vector curvaton model, depending on which region of parameter space is considered, 
T^^'l and r^'^ are certainly very small in vector inflation whereas there is room for them 
to be large in vector curvaton. 

A comparison between scalar and vector contributions can be done by looking at 
the expressions reported in Table 2. Again, one can conclude that the vector contribu- 
tion being dominant or subdominant compared to the scalar one very much depends on 
the parameter space region we want to consider (similarly to the disussion above and 
recalling that the contribution from the scalar field fluctuations in slow-roll inflation is 
at most of order e). 



Table 2. Order of magnitude of ther ratios t^i^/tj^'j^ in different scenarios. 





v.l, S.l 

'''nlI'''nl 


V.2 1 S,l 


V,3 1 S.l 

'''nli'''nl 


V,i 1 s.l 

'''nlI'''nl 


general case 




10-5^3/2,1/2^2 (A) (^n^) ^2^^^^ 






V. inflation 


same as above 


10-5/33/2eV2g2 (A) l^jr^^ 


fie 





V. curvaton 


same as above 


10-5,^3/2,1/2^2 (A) 




r-/33/2,i/2 (-iziP^S 



It is possible to observe that Jnl, as calculated in [55] (see Table 1 therein), and 
T^]^ depend on the same set of parameters (especially if the non-Abelian contribution 
to the bispectrum non-linearity parameter, i.e. /jvl ? is taken into account); this can be 
interesting for different reasons: for example, in the event that the bispectrum and the 
trispectrum were independently known, we could then attempt to restrict the extent of 
the parameter space of the underlying theory; also, it is evident from Table 3 that there 
is enough room in the parameter space of the models we considered for tat^ to be large 
in spite of a small unobservable Jnl- 



Table 3. Order of magnitude of the ratios TjS^^,/ (/jvl ) different scenarios. 





^V.l, (fNA)'^ 
' NL> \-l NL ) 


^V,2 , |'.JVA\2 


^V,3 , (fNA)^ 
' NL' \J NL J 


^V,4 , C.JVA\2 
' NL' \J NL ) 


v.i. 


109£a±« ( H ^2 
\mp ) 


,„e3/2(1+,3) I a\ ( H \^ 







v.c. 


in9 '■^f(l+<3) 


^^r'^e^/^(l + l3) mp m'p 




ir,6 r3j.V2(i + ^) m% jj4 


aip'^ A^ A^ 


/3=/2g2 A^ H A2 




lU g2^5/2 A'i Ai 



In Table 3 f}^£; represents the amplitude of the bispectrum coming from the non- 
Abelian interactions which, as explained in Ref . [55] , can be the dominant contribution. 
It is given by f^f; = 10~^ {Pgc/{1 + P))'^ {m/H)'^ , m ~ mp in vector inflation and 
m ~ Aj ^ in the vector curvaton model. 
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7. Profile of the trispectrum 



The shape of a correlation function such as the bispectrum or the trispectrum is in 
general very peculiar to the specific model or class of models under investigation, this is 
why it is important to be able to carefully analyse these profiles and to find distinctive 
features that can make them easily recognizable. 

As to the bispectrum and for theories where it is isotropic, this task is fairly simple 
given the small number of parameters it depends on (i.e. the external momenta moduli 
ki, /c2 and k^). In theories where it is anisotropic, the bispectrum profile becomes harder 
to analyse: the number of free parameters notably increases with the addition of the 
angles defining the existing preferred spatial directions. We provided an example of 
this in [55], where we proceeded by first analyzing the isotropic part of the bispectrum 
in order to identify the preferred momentum configuration and then, in that specific 
configuration, by studying the effect of the anisotropic part on amplitudes and profiles. 
Producing a plot of the complete shape was not possible though, unless narrowing the 
number of angular variables down, by being more specific about the spatial orientation 
of the Na vector w.r.t. the ki. 

For the trispectrum, the situation is already complicated at the isotropic level. In 
fact, even in this case, the total number of parameters is too large for straighforwardly 
generating a final plot. It is then necessary to restrict the total number of parameters 
with the help of some assumptions about the relative spatial orientation or about the 
moduli of the wave vectors. A list of possible configurations for the tetrahedron made 
up by the four wave vectors, was provided in [34], where the trispectrum is calculated 
for a general Lagrangian in single-field infiation. 

We are going to proceed in a similar fashion and, in analogy with what was done for 
the bispectrum [55] : we analyse the isotropic part first and, in the next section, we plot 
the anisotropic coefficients in a sample momentum configuration. The main purpose is 
to provide a hint of how the existence of preferred spatial directions is responsible for 
modulating the amplitude of the trispectrum, given a particular momentum dependence. 
Our study of the profile by no means intends to be exhaustive, in fact we are 
going to consider only two of the many possible tetrahedron categories (the so called 
"equilateral" and "specialized planar"); also, we will not plot all of the contributions to 
the trispectrum we analysed in this paper, but only focus on the contributions to the 
vector-field trispectrum, arising from vector-exchange and contact-interaction diagrams. 
We leave the task of working out a more systematic and complete shape analysis for 
future work. 

7.1. Equilateral configuration 

It is defined by setting all of the momenta moduli equal to each other [ki = k2 = k^ = 
k^). The variables of the plot can be chosen as x = k^2/^^ V = ^14/^1; which vary 
in the interval [0,2], with k^^/ki = a/4 — — y'^. The plots of the isotropic parts of 
three vector-exchange (from Eq. (58) and its permutations) and of the point-interation 
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v.e.(I) 




Figure 2. Plots of the contact interaction and of the vector-exchange contributions 
in the equilateral configuration. In this and in the next figures, we use the notation 
c.l. = k-^ dx) j^j^j^j^, v.e.{I) = k-^xlx Eti^- = ^"^ x /J x Y^Ls^i 

and v.e.{III) = k^^ x III x Y^Lg^i (the quantities /, //, /// and ti are introduced 
in Eq. (56)). 

contributions to the trispectrum in the equilateral case are provided in Fig. 2. 
The contact-interaction plot is, as expected, fiat on the whole integration domain, since 
it does not depend on kf2 and k^^.This behaviour is very similar to the one of what 
was defined in [34] as Tioc2-, i-e. the local trispetrum term that is proportional to the 
parameter g^^L of local cubic non-linearities. 

The vector-exchange diagrams instead, diverge as k^^ as (k^-fki) — > (limit in which 
the tetrahedron becomes fiat), z = 1,2,3 respectively for (v.e.I), [v.e.II) and {v.e.III). 
This behaviour is similar to what was defined in [34] as T/oci, i-e. the local trispetrum 
term that is proportional ot f%i^. For the vector-exchange plots it is possible to check 
that the amplitudes run from a very large (~ 4 x 10^) and negative number at x, ?/, 2 ^ 1 
respectively for [v.e.I, II, III), to a very large (of the same order) and positive number 
when x,y and z approach 2 (we define z = kf^/ki). This behaviour can be understood 
by looking at the analytic expressions which are characterized, on one hand, by a very 
smooth power law x dependence (that is also responsible for the sign change when 
moving in the interval x = 1 to x = 2), on the other hand by very large numerical 
coefficients. 

We can then conclude that the shapes of point-interaction and vector-exchange diagrams 
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Figure 3. Plots of the contact interaction and of the vector-exchange contributions 
in the specialized planar configuration (plus sign). 



in the equilateral configuration are completely different. 

1.2. Specialized planar configuration 

This is the case defined by ki = = k^^ and by requiring the tetrahedron to become 
flattened. The expressions for kf2 and kf^ are then given by 

^ = \/l + ^±f V^-a:^)(4-y^) (93) 

^ = ^x^ + y'-^T f ^(4 - x2)(4 (94) 

where x = ^2/^1 and y = k^/ki and varying between and 2. The plus and minus signs 
in Eqs. (93) and (94) indicate a double degeneracy in the structure of the quadrangle 
to be expected in the specialized planar case: the plus sign corresponds to a quadrangle 
with internal angles > vr, the minus sign to < tt [34]. The treatment of this case deserves 
some attention in the sense that, while in [34] the plus and minus cases are perfectly 
equivalent given the symmetry of the trispectrum w.r.t. the exchange of any pair of 
its ki {i = 1, ...,4) with each other, this symmetry does not a priori hold in our model 
because of its anisotropy. The two cases will therefore be considered separately in Figs. 3 
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c.i. 




Figure 4. Plots of the contact interaction and of the vector-exchange contributions 
in the speciahzed planar configuration (minus sign). 

and 4. 

Let us begin with the contact interaction diagram. This is independent of fcfg and /cfg, 
it has therefore the same expression for both plus and minus sign cases. The graph 
diverges in the hmit x — > as x~'^ and similarly for y (the expression is symmetric in x 
and y) and it is roughly flat in the remaining part of the domain. 
As to the vector-exchange plots, the situation is very similar. In the plus sign case, the 
diagram (v.e.I) diverg ^ and y ^ respectively as x — j- and y — )■ and is fiat 

for the rest of the domain. The diagram [v.e.II) shows similar divergences in the limit 
a; — )■ and y — 0. As to {v.e.I 1 1), besides for the regions a; — and y — )■ 0, it also 
blows up for X ^ y. This last divergence is due to the fact that (v.elll) goes like k^^ 
and, when x = y, k^-^ = 0. 

In conclusion, for this configuration, we can say that the contact-interaction together 
with the vector-exchange (/) and (//) graphs have similar features in terms of 
divergences and plateau locations. The graph (///) has additional divergences for x ~ y, 
a feature in common with the local trispectra Tiod [34]. 

Let us now come to the last set of graphs: specialized planar configuration with minus 
sign for the vector-exchange diagram (Fig. 4). First of all, we can notice that, while 
for v.e.{I) and v.e.{III) the plus and minus specialized planar cases are different, they 
conincide for v.e.{II), being this independent of k-^2 k-^-^. 
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Figure 5. Plots of the anisotropic part of the trispeetrum from the contribution 
due to vector-exchange diagrams in a sample angular configuration. The analytic 
expression for this plot is cos^ 6 cos^ 5 and was derived from Eq. (97) after normalizing 
to glHtN% [ISO] and setting x* = 1. 



As to v.e.{I), it is singular for x = and y = (the asymptotic behavious is in fact 
given by ~ and respectively as x — )■ and ?/ — )■ 0). The plot shows additional 
singularities in the central region of the plot due to /c^g ~^ 0- to v.e.{II), it has 
singularities for a; — )■ (~ x~^), y — > (~ y~^) and {x,y) — )■ (2,2) (~ {x — 2)~^/^ and 



8. Features and typical level of anisotropy in the trispeetrum 

Let us now consider the complete expression for the trispeetrum, i.e. let us include 
its anisotropy coefficients. Again, we will focus on the non-Abelian trispeetrum 
contributions and, in particular, on the vector-exchange one, given that, in general, 
its amplitude turns out to be several orders of magnitude larger than the one of the 
contact-interaction diagram. In order to discuss an example and provide a plot of the 
anisotropic part of the trispeetrum, it is convenient to pick a well-defined configuration 
reducing the total number of variables to some < 2 number. One possible choice is to 
consider the case where the tetrahedron is fiat and the wave vectors form a rectangle. 
We need to partially fix the orientation of the vectors A^^ w.r.t. the wave vectors. Let 
us pick the following sample configuration 



(y- 2)^3/2). 



N2 




k 



0{i = l 
cos 6, 
cos 9, 



4) 



Ni-k2 







In addition to that, let us assume that all the Na have the same magnitude A^^- In this 
configuration, we have 



= Am = cos^ 9 cos^ 5, Au = 0, (96) 

therefore the the expression in Eq. (58) becomes 

T<^ D g^HtNX [ISO] cos^ 9 cos^ 6 (97) 
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where the expression in brackets includes an isotropic term (which is rotationally 
invariant) 

ISO^I^^y(^IpU + inf:t}j. (98) 

Fig. 5 plots the trispectrum contribution in Eq. (97) normalized to its isotropic part. 
It shows that the trispectrum gets modulated by the preferred directions that break 
statistical isotropy. 

The trispectrum amplitude tj^l was calculated in the equilateral configuration in Section 
6 (without taking into account the angles between the wavevectors and the vectors). 
The amplitude is actually modulated by the anisotropy coefficients, as the previous 
example shows. Notice that in general the bispectrum and trispectrum can be expressed 
as the sum of an isotropic plus an anisotropic parts. This was for instance done in [45] 
for the Jnl parameter in the U{1) case. For our case, the purely isotropic parts of the 
trispectrum are included in the last two lines of Eqs.(30), as far as the Abelian part is 
concerned, and are given by Eqs.(37) and (84), for the non- Abelian terms. The level 
of the anisotropic and isotropic parts can be read from Tables 1 and 2. In particular 
%L ,rNL and Tj^l give the order of magnitude for the level of the corresponding isotropic 
and anisotropic contributions, showing that the two are comparable; on the other hand 
Tj^L quantifies only a pure anisotropic contribution, and it turns out that it can be 
dominant w.r.t. the isotropic part of the full trispectrum. 



9. Trispectrum for /(r) models of gauge interactions 



As anticipated in Section 2, we will now show that it is quite straightforward to extend 
the calculations we performed for / = 1 to cases where / is not a constant. One 
interesting model is the one studied in [35] and also recently discussed in [45] , where the 
field is effectively massless (mo = ^ = 0) so the action (1) for the gauge field becomes 

1 



s' = y c/^xa/^ 



y i^iu pa 



(99) 



/Sf = aAl. The 



where again F^^ = df,B^ - duB^ + g^e" 

Let us introduce the fields A" and A", related by the equations A1 
A" are the physical fields. 

We can expand the perturbations of in terms of creation and annihilation operators 
in the usual way 



<5i«(r?,f) = /-^ [e^{q)af5Al{r^,q) + h.c. 



(100) 



li f = /oct", with a equal either to 1 or —2 (/o is a constant), it is possible to prove 
[45] that the equation of motion for 6A'l is the same as the one for 6B^, where by 6B^ 
we mean the transverse mode function in Eq. (26). This is equivalent to saying that, 
under the assumption a = 1,-2, the physical gauge fields are governed by the same 
equation of motion as a light scalar field in a de Sitter space and so they generate a 
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Figure 6. Plots of the isotropic functions of some of the vector-exchange contributions 
in the equilateral configuration, for the / ~ a^^ model. In this and in the next figures, 
"v. 6.(1, II, III)" represent the isotropic functions associated with the very last term in 
square brackets in Eq. (Ill); "v.e.(new)" represents the isotropic function associated 
with the A:ii44, ^2244, ^1133 and /c2233 terms in the second line of Eq. (111). 



scale invariant power spectrum. 

We are going to sketch the calculation of the trispectrum in this theory. 

The general expressions in Eq. (30) obviously still hold as well as the compact structure 

of the Abelian terms, Eq. (34), except that the power spectrum in Eq. (16) reduces to 

Pti = T^'P+. (101) 
having gauged the longitudinal modes away. 

Let us now have a look at the non- Abelian part and see how Eqs. (36), (46) and (62) 
change when we switch to massless /(r) Lagrangians. First of all, we need to set 
n{x*) = 0. The interaction Hamiltonian to third and fourth order are the same as in 
Eqs. (39) and (40), except for an extra P{t) factor. The anisotropy coefficients Ieee, 
"Pijki^^^ and T^^^^ that survive after setting the longitudinal mode to zero do not 
change. On the other hand, the wavefunctions for the gauge fields are now given by 
6B = 6A/f = 5i?^//oa". The new trispectrum therefore differs because of extra scale 
factors inside and outside the time integrals, which in general imply a different power 
of in the final results and a different momentum dependence in the isotropic part of 
the expressions. 

The three non-Abelian vector contributions in the / = 1 case can be schematically 
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Figure 7. Plots of the isotropic functions of some of the vector-exchange contributions 
in tlie specialized planar configuration (plus sign), for the / ~ model. 




Figure 8. Plots of the isotropic functions of some of the vector-exchange contributions 
in the specialized planar configuration (minus sign), for the / ~ model. 
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written as follows 

{5B%^,2 ^ {SB') J drj {SBf = {SB^f J dry (SB^f , (102) 
(55^),.,. ~ {6B)'J dr]'{6B)'J dr]"{6Bf 

'5B^) ' [ dv' (55^) ' / rfr/" (5B^) \ (103) 



(55^),.,. ^ {6B)^J dr^{5Bf = {SB^yj dr^(5B^)\ (104) 

respectively from Eqs. (36), (46) and (62), where {5B^) = {5B^(ki)5B](k2)5Bl{h)6Bf(ki)) 
and we have omitted all the gauge and vector indices, as well as complex conjugate sym- 
bols. Let us see now how Eqs. (102)-(104) change if / = f^a^ (a = 1, —2) 

{5B%^,, ^ {5Bf j dvf m' ^ (^^) / dva'- (^^) , (105) 

{5BXe. ^ {5Bf j drlf{5Bf j dri'p{5Bf 



f^^V (^ I. (106) 



'6B^V r . o„, /6B^^^ 



{6B%,,, ^ dvf ^ (^^j / dr^a- (^^j . (107) 

Using a = {—Hri)~^ in the previous equations, we get 



{SB')une2^Ht''{SB^i-vTy I dv{SB^fi-vr, (108) 

X J dr]" {6B^y {-r]"r, (109) 

(SB^),,,, ~ i/f (55^(-r/*)") Y (^5'^)' (-^)'"- (110) 

Let us now consider more in details the a = —2 case for contact-interaction 
and vector-exchange contributions. The expressions for the anisotropy coefficients are 
respectively given by 



J.EEEEE ^ kihi{ki ■ ks-ki- kf^h ■ kf2)[SijSki - Sijlck^kiA - Sijhsks - 5kiki2kj2 - Skikikji 
+ ^ijkkskuks ■ k^ + Skikiikj2ki ■ /c2 + ku4^4 + A;2244 + ^1133 + ^2233 ~ ^2234^3 " ^4 
- ^1134^:3 ■ h - kuMh ■ h + /^i233^i ■ h + ki2Mki ■ hh ■ h] (HI) 

and by Eq. (63), for one of the possible permutations. These expressions are more 
comphcated w.r.t T/jj." in Eq. (53) and T^^^ in Eq. (78) for the longitudinal modes. As a 
result, when studying the shape of the trispectrum, for the isotropic functions appearing 
in it, several diagrams need to be taken into account, one for each term in T^^^i^^ 
Tl^ki^^^ ■ For comparison with the / = 1 case, we plotted the isotropic functions 



29 



associated with the very last term in square brackets in Eq. (Ill) (see "— w.e(/)", 
"— w.e(//)" and "— w.e(//J)" in Figs. 6, 7 and 8). By comparing these plots with the ones 
in Figs. 2, 3 and 4, it is evident that they have very similar shapes. On the other hand, 
when we consider the isotropic functions associated with terms that are not present in 
the / = 1 case, several differences arise in the plots; we provide a sample in Figs. 6, 7 
and 8 with the ^^v.e.{new)" plots, which represent isotropic functions associated with the 
^1144, ^2244, ^1133 and A;2233 tcrms in the second line of Eq. (111). We verified that similar 
observations can be made concerning the shapes of the contact-interaction contributions. 

10. Overview and conclusions 

We have calculated the trispectrum of curvature perturbation in an SU{2) model of 
gauge bosons coupled to gravity during inflation, where the latter is driven by a slowly- 
rolling scalar field, Eq.(l). The SN formalism has been employed to relate the curvature 
perturbations to the perturbations of all the primordial fields. 

The trispectrum is made up of contributions due both to the scalar, Eq. (7), and to 
the gauge fields, Eq. (30). The latter can be distinguished in "Abelian", i.e. retrievable 
in the Abelian limit of the theory, Eq. (34), and "non- Abelian" , i.e. specific to our 
model, Eqs. (36), (58) and (62). It is worth to point out that as a by-product of our 
results, the trispectrum in the pure Abelian case has also been computed for the first 
time. All of the vector field contributions can be written as the product of an isotropic 
function of the momenta moduli (fci, k2, k^, k^, k^2 and k^^ and of the time labeling the 
initial hypersurface 77*, times anisotropy coefficients that depend on the angles defin- 
ing the direction of the gauge fields w.r.t. each other and to the wave vectors. The 
non-Abelian contribution, calculated using the Schwinger-Keldysh formalism, is made 
up of a term that is proportional to the bispectrum of the gauge fields, Eqs. (36), and 
terms that originate from the trispectrum of the gauge fields, Eqs. (58) and (62). The 
latter can be diagrammatically represented by two categories, depending on the spe- 
cific interaction Hamiltonian term involved: vector-exchange and contact-interaction. 
They present a similar analytic expression in terms of powers of the SU{2) coupling 
constant and of the Hubble parameter H^,, they are though very different in terms of 
their anisotropy coefficients and of their momentum dependence. 

We have verified that, in the case where the longitudinal and the trasverse com- 
ponents of the wave vectors have exactly the same evolution and the model does not 
violate parity, the trispectrum becomes isotropic in all its contributions, except for the 
vector-exchange one. We also showed that the non-Abelian part of the trispectrum van- 
ishes in the case where all the gauge fields are aligned along a unique spatial direction. 

We have studied the amplitude ttvl of the trispectrum. We neglected all the gauge 
and vector indices (Table 1), in order to provide an approximate evaluation of tnl in 
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terms of the parameters of the theory. We showed that the impact of these parameters 
on suppressing or enhancing the amphtude of the trispectrum very much depends on the 
specific subset we decide to pick in their parameter space (it depends also on the choice 
of a specific gauge field configuration, i.e. on the magnitude of their initial value and 
their spatial orientation). We also noticed that, when comparing the contact-interaction 
and the vector-exchange contributions (which present a similar parametric dependence 
from the SU (2) coupling and from the Hubble parameter -ff*), the latter are generally 
several order of magnitude larger than the former. This is of course true unless we are 
in a configuration for which the anisotropy coefficients somehow reverse this situation 
ending up by suppressing the vector-exchange contribution. 

We finally carried out some shape analysis, focusing on the contribution due to the 
trispectrum of the gauge fields. As specified above, each contribution to the trispectrum 
is written as the product of a term that is invariant under rotations of the tetrahedron 
times some anisotropy coefficients. The study of the shape was performed in two steps: 
first the rotationally invariant part was analysed for two distinct momentum configura- 
tions, the equilateral and the so called specialized planar configurations. Figs. 2, 3 and 
4, then the anisotropic coefficients were evaluated in a sample angular configurations. 
Fig. 5. For the equilateral case we found many similarities between all our plots and the 
shape of a local trispectrum, whereas, for the specialized planar case, our plots turn out 
to be easily distinguishable from the local plots except for one of the vector-exchange 
graphs {v.e.III). Also, the plots for point-interaction and vector-exchange have some 
common features in the specialized planar configuration whereas they turn out to be 
completely different in the equilateral case. In Fig. 5 we showed how the amplitude of 
the trispectrum is modulated by the presence of the anisotropy coefficients. We also 
verified that the amplitudes of the anisotropic and isotropic parts of the trispectrum can 
be of the same order of magnitude. This last analysis was not meant to be exhaustive 
but merely provide the reader with an idea of how the shapes can be studied and how 
anisotropy can leave a strong imprint on them. 

Studying the generalization of our calculations to other gauge symmetries is an 
interesting possibility and it would translate into different results w.r.t. the SU{2) case. 
In general, the expression of the interaction Hamiltonian would change and therefore 
the bispectrum and the trispectrum of the vector bosons are expected to be different. 
Consider for example a generic SU{N) group, with N > 2; one realizes that, since the 
number of vector bosons is equal to the number of generators of the group, the number 
of preferred spatial directions would increase as we go to higher A^. This is an important 
feature that allows to distinguish among different symmetry groups. 
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Appendix A. Anisotropy coefficients for the trispectrum mixed 
contributions 

We list the coefficients appearing in Eq. (29) 



h ■ Na) {h ■ Nc) k-h + {k ■ No) {h ■ Nc) k ■ h) + (l ^ 2) + (3 ^ 2) 
(2 [Na ■ Nc) {k ■ N^b) {k ■ A')) + (2 ^ 1) + (2 ^ 3)" 
[h ■ A- [2 ■ Na) {h ■ N^b) {h ■Nc)+2 {h ■ Na) {h ■ N^b) {h ■ N,) 

k ■ Na) {k ■ N^b) {k ■ iVe) k-h- {k ■ Na) {k ■ N^b) {k ■ N,) k ■ h 

+ (2 ^ 1) + (3 ^ 2)]], 

k ■ Na) {k ■ N^b) {k ■ N,) [k ■ k) {h ■ 



+ 



(A.l) 



AS) 
'■III 



a'ab^a'ce 



e e 



k ■ Na) {k ■ N^b) {k ■ Nc) {k ■ k) {k ■ A')) + {l ^ 3) + {2 ^ 3)] , (A.2) 



r(8) 
'llE 



r(8) 
'■EEl 



^a'ab^a'ce 



^a'ab^a'ce^^^^^ . ^e^ ^ ^J^^ . ^j^^ . ^ ^J^^ . . ^ . 

^(2 ■ A-) [k ■ Na) {k ■ N^b) {k ■ iVe) ) + (1 ^ 2) 

( ((fci ■ Na) {k ■ iVc) + {k ■ Na) {k ■ N,)) [k ■ N^b) {k ■ A') k ■ k 

+ (1 ^ 3) + (2 ^ 3)]], 

4 [N^b ■ A') {k ■ Na) {k ■ Nc) 

k ■ N^b) {k ■ A-) {{k ■ Na) {k ■ iVe) + (^3 ■ Na) {k " N,)) k " k 

+ (2 ^ 1) + (2 ^ 3)" 

(2 [k ■ A') [k ■ Na) {k ■ N^b) {k ■ N,)) + (1 ^ 2) + (2 ^ 3) + (1 ^ 3) 
[[N^b ■ A') k ■ k {{k ■ Na) {k ■ N,) + [k ■ Na) {k ■ N,))) + (1^2) 
[Na ■ N^b) {k ■ N,) [k ■ A') + {Nc ■ N^b) {k ■ Na) {k ■ A 



(A.3) 



+ 



+ 



(A.4) 



where N^b^ {Nib, N^b, Nib). 

As to the coefficients lajs-y appearing in Eq. (36), they can be easily derived by replacing 
Nb to N^b in the expressions of where we defined A* = M^N^. 
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Appendix B. More details on computing the analytic expressions of 
vector- exchange diagrams 

We report the expressions of the functions A, B...P introduced in Eq. (47) 



A = i^-lGk'^ + kik2X*'^) cos 

B = A[ki^ ks, k2 k^] , 

C =4:k {ki + /C2) X* cos 

D = C[ki ^ ks, k2 ki 



{ki + k2) X* 



4k 

{ki + k2)x*' 
ik 



— Ak {ki + /C2) X* sin 



E = l8k (kf2 + k^ + ki) — k^-^k^k^x* j cos 



+ (-16A;2 + kik2X*^) sin 

(/Cf2 + ^3 + h) X* 



(fci + k2)x* 
Ak 

(fci + k2) X* 
4k 



4k 



+ 2k {k{2 + ^3 + ^if ^* sin 
F =E[k^^ ki,ki k2], 
G = 2k {k{2 + ^1 + ^2)^ X* cos 



(/Cl2 + ^3 + ki) X* 
4k 

"(A;f2 + ki + k2)x* 



4k 



+ i-Sk + ^1 + ^2) - ^12^1^22^ 



*2 



sm 



{k{2 + + k2) X* 
4k 



k^ ~\~ k-^2 ^-|22^2 ~t~ ^i"2^2 ^~ ^2(^12 ~^ ^2) ~t~ ^l(^]"2 ~^ ^2)) *^ 



2\\ *2 • 

2) x a 



(B.l) 
(B.2) 
(B.3) 
(B.4) 

(B.5) 
(B.6) 



(B.7) 

{k-i^2 + + k2) X*' 
4k J ' 

{k^2 + h + ^2) X*' 



4k 



M = (/c? + k% + A;^~2^2 + A;i2^2 + ^2(^12 + ^2) + kiik}^ + k. 

N = M[k3^ ki,k4 k2] , 
P = L[ki ^ k-s, k2 -> ki] . 

The anisotropy coefficients introduced in Eq. (56) have the following expressions 



(B.8) 
(B.9) 



ti = kiks (Jvi ■ k 
t2 = kiki (ki ■ k 



"12 

ts = k2k3 {k2 ■ ki2 
ti = k2ki (ji2 ' ^12 
= kik2 (J^i ■ A;f3 
= kiki (ki ■ kf^ 



tl = ^2^3 (^3 ■ ^13 

ts = kski (k3 ■ kf^ 
tg = kik2 (ki ■ 



ki ■ k 
ki ■ k 
ki ■ k 
{ki ■ k 
{h ■ k 
{ki ■ k 
{ki ■ k 
{h ■ k 
[h ■ k 



h ■ h 
h ■ h 
{h ■ h 
{k2 ■ ki 
{k2 ■ ki 
{k2 ■ ki 
(Jt2 ■ ki 
ko ■ k^ 



(jti ■ ky^ 



h ■ ^12 



(Jvi ■ k 

(^2 ■ k^'r 

(Jvi ■ k 
{k2-k 



ko, ■ k 



14 
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(B.IO) 

(B.ll) 
(B.12) 
(B.13) 
(B.14) 
(B.15) 
(B.16) 
(B.17) 
(B.18) 



tio = hh (ki ■ kf^'j [ki ■ kij (k2 ■ h) [h ■ kf^'j 
tn = kiks (ki ■ ^12) (^1 ■ ^2) {h ■ ^4) (ks ■ k^^) 
ti2 = kiks (ki ■ ^12) (^1 ■ ^2) (h ■ h) (h ■ ^12 



(B.19) 
(B.20) 
(B.21) 



Let us now list all the scalar products appearing in the equations above 



A; -|- A; -|- fc, 
12 ^ 2fcifc,-„ 



ki ■ k 



13 



ki ■ kf^ 



14 



k2 ■ kf. 



k « ~\~ k-i k-i 

12 1 1 



12 



2k2k 



k?. ■ k 



13 



^k It ~\~ k A k k ~\~ k i-. 



k/i ■ k 



14 



fc4 -kf- 

14 1 



2fc4fc • 



k?. ■ k 



12 



h'^ |_L2 _9L2_i,2_i%2 
12 14 2 1 3 



ko ■ k 



13 



1,2 i_i,2 _ 1.2 _ 1.2 n 1.2 

14 12 1 3 4 



i-2 i_i,2 _JU2_9i,2_i,2 
1'2^ 1'4 1 2 3 
2fcifc3 



/c2 ■ kA 



h.2 \U2 _JU2_t,2_9i,2 
14 12 1 3 4 



2/^2 /i;4 



~T i,2_r.2 _i.2 

7 i, . _ ^^3 '^14 ^=2 
'*^2 «;i4 - 2fc2fcj^ 



/c/i ■ k 



12 



'2ik^ k-^^ 



ki ■ k2 



^l+^2~^12 
2k^ko 



ki ■ ks 



k A ~\~ k — k ' — k ' - 
4 ' 2 12 14 



2A:ifc.3 



ki ■ k^ 



1.2 _h.2_h.Z 
'^14 "^1 '^4 

2fci kA 



k-i ■ kd 



1.2 _h.2 _h.2 
"-12 "^4 '^3 

2fc4fc3 



kp. ■ k 



14 



h. 2 _h. 2 _U2 
2 3 14 

2^3 ^-] 4 



A;2 ■ k-?. 



UZ _iu2_iu2 
14 2 3 

2fc2A:3 



where ki = \ki\, ki = ki/ki, k^- = ki + kj, kij = \ki + kj\, i and j running over the four 
external wave vectors. 



Appendix C. Complete expressions for functions appearing in 
point-interation diagrams 

We provide here the expressions for the coefficients appearing in Eq. (79) 

Qeeee =x*^[- k{ek{ + k^^kl -kl + ekl + k^k^ki - ?,k^klk^ki + k^klki 

- kl + k^kskl + k^kl + k^k2{kl - kk^k^ - Sklk^ - Skskl + kl 

+ k\k3 + k^)) + kl{ - 1 + k\k-i + ki)) - Sk'^klikski + k2{ks + ki)) 
+ kl{-l + k\k2 + A;3 + k^)) + k^ki{kl + A;| - ?,klki - Sk^kj + kf 

- 3kl{k3 + ki) + P{k2 + k3 + ki) - 3k2ikl + ^k^ki + kl) - k^k^ki 

+ kiiks + ki))))] + x*^[k'^{kfk2k3ki + kik^kiikl + kl - 8/^2^:3/^4 + kfj 
+ k^ikski + k2{k3 + ki)) - 3kl{klkl + k2k3ki{k3 + ki) + kl{kl + kski 
+ kl)) + ki{klk3ki + kl{k3 + ki) - 3klk3ki{k3 + ki) + k3ki{kl + kl) 

+ A;2(A;3 + klki - 3klkl + kl + kski^k^ + kl))))] - 3x*'^ k^k^klkl, (C.l) 

Aeeee = k{kl + kl + A;3' + kfj x*^ (C.2) 

Beeee = k^ -k'^ [kski + k2 [k3 + ki + ki (/c2 + ^3 + ki) x*^ + kik2k3kiX*^)) , (C.3) 

Ceeee = kx* {k^ - {k2k3ki + kl {kski + /c2 (/C3 + h))) x*^) , (C.4) 

Deeee = - kx* {k^ - {k2k3ki + ki {k3ki + A;2 {k3ki))) x*^) , (C.5) 

Eeeee = k^-k'^ {k3ki + k2 {k3 + ki) + ki {k2 + k3 + ki)) x*^ + kik2k3kiX*^, (C.6) 

Feeee =k(kl + kl + kl + kt) x*\ (C.7) 
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